Unusual decoherence in qubit measurements with a Bose-Einstein condensate 
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We consider an electrostatic qubit located near a Bose-Einstein condensate (BEC) of noninter- 
acting bosons in a double-well potential, which is used for qubit measurements. Tracing out the 
BEC variables we obtain a simple analytical expression for the qubit's density-matrix. The qubit's 
evolution exhibits a slow (oc 1/Vf) damping of the qubit's coherence term, which however turns to 
be a Gaussian one in the case of static qubit. This stays in contrast to the exponential damping pro- 
duced by most classical detectors. The decoherence is, in general, incomplete and strongly depends 
on the initial state of the qubit. 



PACS numbers: 03.75.Gg, 03.65.Ta, 03.65.Xp 



I. INTRODUCTION 

Recent progress in quantum information technology 
has lead to significant technological and theoretical ad- 
vances in measuring and controlling the state of a two- 
level quantum system (qubit). Devices used for this 
purpose include point-contact detectors and single elec- 
tron transistors [l|, 0, Q where the magnitude of elec- 
tron current is used to determine the qubit's state. Re- 
cently, more sophisticated hybrid systems which combine 
a charged qubit with microwave resonators or ensembles 
cold polar molecules were proposed 4, 5]. In addition to 
technological benefits such hybrids offer an insight into 
fundamental physical phenomena, such as decoherence. 
The decoherence is present for any microscopic system 
(e.g., a qubit) interacting with a macroscopic device, 
characterized by a large number of degrees of freedom 
and a dense distribution of energy levels. As a result, 
an initial state of a qubit is expected to be rapidly (ex- 
ponentially in time) converted into a statistical mixture, 
so that the information stored in the qubit is erased Q. 
For example, this has been explicitly demonstrated for 
qubit measurements with a point-contact detector shown 
in FigQJ,, where a macroscopic current flowing into the 
right reservoir across the potential barrier is modulated 
by the qubit's electron |7|]. 

In this letter we study measurements in a hybrid sys- 
tem, consisting of an electrostatic qubit placed in closed 
proximity to a non-interacting BEC trapped in a sym- 
metric double-well potential. The qubit is represented 
by an electron in coupled quantum dots (FigfTJa), while 
confinement of the BEC can be realized, for example, 
by means of a quasi-electrostatic optical dipole trap pro- 
duced by two crossed laser beams @ . Since the trapping 
occurs due the interaction of the induced atomic dipole 
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FIG. 1: (Color on line) a) Point-contact in which electrons, 
initially in the left reservoir, tunnel across a potential barrier 
modulated by the presence of an electron in one of the two 
coupled quantum dots (qubit); b) Atoms of a BE condensate 
tunnel across a barrier of a symmetric double-well potential 
modulated by a qubit. 



moment of neutral atoms and the far-detuned optical 
field [3], additional electric filed induced by the electron 
would change the barrier height. This, in turn, would 
modulate the atomic current just as the presence of an 
electron in one of the dots modifies the current of the 
point-contact detector shown in FigQJ,. There is, how- 
ever, an important difference as only a single level (zero- 
width band) is available for the tunneling atoms, which 
raises the question of what type of decoherence, if any, 
would experience the measured qubit? 

With the number of carriers macroscopically large, but 
only one level existing in each of the reservoirs, the ques- 
tion cannot be answered without a detailed analysis. On 
one hand, the bosons are moving independently, and one 
could expect their effect to be similar to that of a single 
boson which, as is easy to show, does not produce deco- 
herence. On the other hand, it is not clear whether the 
large number of uncorrelated degrees of freedom in the 
detector will not have an averaging effect on the qubit 
thus causing the off-diagonal elements of its density ma- 
trix to disappear. 



In the following we will show that the density matrix 
of a qubit coupled to a BEC will undergo an evolution 
which is not be described by Bloch-like equations similar 
to those arising in the case of a point-contact detector [7| . 
However, the BEC model shown in FigflJ) can be solved 
exactly. We will demonstrate that the rate of decoherence 
is extremely (non-exponentially) slow and, unlike in the 
case of the point-contact, its amount strongly depends 
on the choice of the qubit's initial state. 



is the state corresponding the qubit localized in one 
of the quantum dots, q = 1,2, and n = 0,1,2,... is 
a number of bosons contained in the right well. For 
the corresponding probability amplitude ipqn (t) we write 
4>qn(t) = (q\ (4>n I sxp(—Ht)\qo)\(f>o) state from which the 
reduced density matrix of the qubit is obtained by tracing 
out the BEC states, 



°OT'(*) =^2i>qn(t)^q> n {t) . 



(4) 



II. DESCRIPTION OF THE MODEL 

Consider an electrostatic qubit interacting with a BEC 
consisting of N atoms initially trapped in the left well of 
a symmetric the double- well structure (Fig. [TJd). The 
entire system can be described by the tunneling Hamil- 
tonian H = H c + H q + H mt , where the three terms corre- 
spond to the condensate, the qubit and the interaction, 
respectively, and 

H c = E (c{c L + c R c R ) - fl{c[c R + c R c L ) (la) 

H q = E x d[di + E 2 d\d 2 - Lu (d\d 2 + d\di) (lb) 

H int = 5Vt d\d 2 {c\c R + c R c L ). (lc) 

Here c L R is a boson creation operator in the left (right) 

reservoir, and d\ 2 is the fermionic creation operator for 
the qubit. Tunneling between the reservoirs is suppressed 
when the electron is in the dot nearest to the barrier, so 
that 50, = O - tt' > 0, Fig|Tb. 

For a static qubit trapped in one of the quantum dots, 
(wo = 0), each boson of the condensate oscillates between 
the reservoirs with the Rabi frequency O or fi'. Thus, the 
probability of finding n bosons in the right reservoir at 
the time t is given by 



P n (t) = ( N j cos 2{N - n \nt) sm 2n (nt), 



(2) 



if the qubit's electron occupies the level E\ of the nearer 
dot (or by the same expression with Q — » Q' when it 
occupies the level E 2 ). In the interesting case when 
the tunneling rate for each atom is small but the num- 
ber of atoms is large, we put N — ► oo, while y/~NQ — > 
K =const, thus maintaining a finite current into the right 
reservoir. For small times, t — At << f2 _1 Eq. @ 

yields P„(At) « (KAi) 2n e"' KAt ) 2 /n! and, in particular, 
Po(At) rs 1 — (nAt) 2 . Such a non-Markovian behavior 
of the BEC is in contrast with large fermion reservoirs 
(Fig. [T]a,), where one finds 1 — Po(Ai) oc (At), which is 
typical for a Markovian process. 

Consider now the behavior of a dynamic qubit, ujq ^= 0, 
subjected to a measurement with such a non-Markovian 
(BEC) detector. The wave function of the entire system 
can be written as \^>(t)) = Y, q , n VvWItf)!^") where 



I iii \ i+in\ \ C l) \ c R/ \v)cond 

\q)\(t> n ) = d\\Q) qub 

\/\JS — n)\n\ 



(3) 



Putting, for convenience, Eq = one easily finds from 
Eqs. {1} that H q + H lnt commutes with H c and the evo- 
lution operator U(t) = exp(— Ht) can be factorized, 

U(t) = e -^t e ^(H,+H lat )t s C/c ( t )[/ q . nt(t ). ( 5 ) 

Operators U c (t) and U q i n t(t) can be written in a simple 
form in basis of the eigenstates of the BEC Hamiltonian, 



\4>n) = 



(4+4r-"(4-4no) 



y/2 N {N-n)\nl 
(N -2n)n\4, n ). 



(0) 



(q, (f) n | U c (t) | q',(l>n>) = exp[-i(N - 2ri)Q,t]5 nn <8q q . 



such that H c \4>. 

also 

(g|(^|C/ q i„t(t)| g ')|^) = (q\U(t,s n )\q')5. 



Indeed, we have 
and 



(7) 



where e n = (2n — N)6fl. It is readily seen that U(t, e n ) = 
exp[— i(H q + e n d 2 d 2 )t] represents the evolution operator 
of an isolated asymmetric qubit with the level displace- 
ment e' n = e n + E 2 — Ei (asymmetry parameter), whose 
matrix elements are easily found to be 



cos (cut) — i— sin (u>t) 

2to 

e' 
cos (u)t) + i— sin (tot) 
2u> 



e -fe'„*/2 



-i— sin (tot) e 



—ie' n t/2 



(l\U(t,e n )\l) 
(2\U(t,e n )\2) 

{l\U(t,e n )\2) = :■■ 

where w = w{e n ) = ^J{e' n /2) 2 + Wq is the qubit's Rabi 
frequency and (2|l7(t,e n )|l) = (l|L>(i,e„)|2). 

For the reduced density matrix of the qubit in Eq. IJ3} , 
with the help of Eqs. (gj)-© we find 



(8) 



<?qq' (t) 
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(9) 



where 

' qq' v^, Er, 



cr„ 



(g|C/(t,£ n )|go)(9o|f>" 1 (*,£n)k')- (10) 



is the density matrix corresponding to the unitary evo- 
lution of an isolated asymmetric qubit. With the initial 
state of the BEC given by \(f> ) = (l/\/M)(c[) Ar |0) we 
then find 



N\ 



2 N n\{N-n)\ y^N 



(N -2n) z 

e -^^- (ii) 



w here we have used the Sterling formula K\ ~ 
\/27rKK K exp(— K) to evaluate the factorials. 

Now we assume that the qubit's coupling with each in- 
dividual atom of the condensate (SCI) is weak, but its in- 
teraction with the entire condensate is considerable, and 
so is the variation of the BEC current (ex y/N Sfl), in- 
duced by the qubit. Then taking the limit 

2V-*oo, with VnSQ -» a = const (12) 

we replace the sum over n in Q by an integral, J^ n — > 
fcfe/(2£Q). This yields 



/OO 
a qq ,(t,e)exp(~e 2 /2a 2 ) 
-oo 



rf; 



27ra 



(13) 



where a qq > (t, e) = [a qq ' (t, e) + o qq * (t, — e)]/2. 

In the following we will consider only the case of a sym- 
metric qubit, E\ = E 2 10]. Then for an initial qubit's 
state \qo) = (ad\ + bd 2 )\0) , we obtain from Eq. (fit)]) 



cr u (t,e) = \a\ 2 + (\b\ 2 - \a\ 2 )[l - co S (2Lot)}(a J 2 /2tu 2 ) 
- lm(ab*) s'm(2ujt)(uj /uj) , (14a) 

a 12 (t, s) = i{\a\ 2 - \b\ 2 ) sm(2Lut)(aj Q /2uj) + ab* cos(2ujt) 
+ Re(o6*)[l - cos(2Lut)](uj /uj) 2 , (14b) 



with u> = w(e) = \J(e/2) 2 + w 2 ) and <7 2 2(£) = 1 - 0ii(*), 
5-21 (*) = 5f 2 (*). 



III. DECOHERENCE OF QUBIT DUE TO ITS 
INTERACTION WITH THE BEC 



The simple form of Eqs. (fT3")) - ([i"4"l) allows for an easy 
analysis of limiting cases. Indeed, the strength of inter- 
action with the BEC, a, enters Eq. (fT3f only via the Gaus- 
sian cut-off factor exp(— e 2 /2a 2 ). The factor determines 
the number of asymmetric configurations contributing of 
the qubit's evolution and, therefore the perturbation in- 
curred upon the qubit by the BEC. (Note that when the 
interaction vanishes, a — > 0, the Gaussian becomes nar- 
row, and we recover the unperturbed evolution of the 
isolated qubit.) 

It is readily seen that in the large time limit t — > 00 
the contributions for the rapidly oscillating terms in 
Eqs. (ri"4"a V (|14b|) vanish. Evaluating the remaining inte- 
grals analytically shows that as t — > 00 the density matrix 
of a qubit tends to a steady state a st given by 



(V^/2)z exp (z 2 ) erfc (z) (\b\ 2 - \a\ 2 ) , 



'12 



^n zexp (z 2 ) erfc (z)Re(ab*) 



(15) 



where z = \/2uJo/a and erfc (z) = -j= J^° exp(—t 2 )dt is 
the complementary error function. 

To evaluate the speed with which this steady state is 
attained we note that at large t the phase of the sines 
and cosines in Eqs. (fT4")l develops a stationary region of 



the width Ae = (Aujo/t) 1 ^ 2 centered at e = 0. Once 
Ae becomes small compared to the width of the Gaus- 
sian in Eq. (|13|) . i.e. for t >> ojq/o 2 , the contribution 
from the stationary region becomes proportional to Ae 
causing the time-dependent part of o qq ' in Eqs. (|14[) to de- 
cay as 1/y/t. (For a recent discussion of non-exponential 
decoherence behavior expected in other systems see, for 
example, Ref.|13. fl3l|.) Explicitly, for cj ^ 0, the sta- 
tionary phase method yields: 



au(t) ^ a s A + ] ^ [(\a\ 2 - \b\ 2 ) cos fat + \) 



2Im(a&*)sin ( 2uo a t + - 



(16a) 



oia (i) ~ erf* + ^ s ^- [(\a\ 2 - \b\ 2 ) sin fat + \ 



+ 2Im(a6*)sin(2w t + 



(16b) 



where a st is given be Eq. (fTS"]) . Figure [2] demonstrates 
that Eq. (fl6|) (dot-dashed curve) coincides to graphical 
accuracy with the exact result (fTUJ) (solid curve) except 
at very short times. 

Equations (fTS")) . (|16[) which describe the qubit's deco- 
herence generated by the BEC employed as a measure- 
ment device represent our main result. The qubit's be- 
havior is very different from that of a qubit interact- 
ing with electronic reservoirs [7|, Fig. [TJ or in a general 
with any Markovian environment, whose effect can be 
described by Bloch-like equations [lj, [lfl [lg|. Indeed, 
it follows from Eqs. (TIT)]) that the relaxation to the fi- 
nal steady state is extremely slow, obeying the power 
law oc \j\ft. One exception from this rule is a static 
qubit (loq — 0) for which the stationary region van- 
ishes so that from Eqs. fTS")) . (I14p . one easily obtains 
o~i2(t) = ab* exp(— a 2 £ 2 /2). In contrast, in a Markovian 
environment, a static or dynamic qubit undergoes an ex- 
ponential relaxation to the final statistical mixture. 

It also follows from Eq. (|15| that , in general, the 
qubit's decoherence in the steady state is incomplete and 
its density matrix is not converted into a statistical mix- 
ture, a st ^ diag (1/2, 1/2) as would be the case for a 
point-contact detector. Rather, complete decoherence is 
achieved only in the weak coupling limit (a — ► 0) [ll| 
and only for the initial conditions corresponding to Re 
(ab*) = 0. For a weak coupling, the dependence on 
the qubit's initial state can be understood in a follow- 
ing way. A real part of the qubit's off-diagonal density- 
matrix element Re<7i2(£) can be written as Heai2(t) = 
(1/2) (*(t)|<7+ 1*(£)), where q+ = d\d 2 + d\d x . If the 
qubit's levels are aligned (E\ = E 2 ), the operator q + 
commutes with the total Hamiltonian, Eq. ([1]), in the 
limit of 6fl — ► 0. As a result Recri2(£) — Re 1712(0). 
Therefore, the subspace of the qubit's states, correspond- 
ing to Im (112(0) = is effectively decoherence free [17| . 

In strong coupling limit, a — > 00, in Eqs. ([T4"k )- l|14b[) 
we only need to retain the terms which do not vanish 
for |e| — > 00, thus compensating for the factor a -1 in 
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FIG. 2: (Color on line) a) The qubit's density matrix element 
an vs. time t for |g ) = [(l + i)|l) +i\2)]/V3, a = 1, iV -> oo, 
Eq.JTSJ, (solid), A = 100, Eq.©, (dashed) and the stationary 
phase approximation, Eqs. (|16p (dash-dot). Horizontal line 
shows the large time asymptote (|15[1 ; b) same as (a) but for 

0"12- 



Eq . (fT3")) . Accordingly, the off-diagonal density-matrix el- 
ement would disappear at all times for all initial qubit's 
states, <7i2(£) — > 0. However, the result is not the sta- 
tistical mixture, as in the case of weak coupling, but 
<j(t) = diag(|a| 2 , \b\ 2 ). This corresponds to the so-called 
pure dephasing for a static (luq = 0) qubit [18[ whose 
diagonal density-matrix elements remain constant while 
the off-diagonal elements vanish. 

Finally, the finite size effects for a condensate with a 
large but finite number of atoms are shown FigO These 
manifest themselves as an onset of irregular oscillations of 
the qubit's density matrix (dashed lines in Fig|2|). The 
oscillations appear at times comparable with the Rabi 
period of an individual atom in the double well poten- 
tial, t w 27r<5f2, prior to which the qubit's evolution 
agrees with that in the presence of an infinite condensate. 
Mathematically, the effect occurs when the period of ever 
faster oscillating terms in Eqs. (fT4b )- (|14b[) becomes com- 
parable with the separation between the energy levels of 
the condensate and Eq.© ceases to be a valid Riemann 
sum for the integral (fT5|) . Physically, the qubit begins 
to be affected by the size of the condensate at times of 



the order the Poincaire recurrence time of the latter, i.e., 
when the escape of the atoms into the right well can no 
longer be considered irreversible. 



IV. SUMMARY 

In summary, we have demonstrated that continuous 
monitoring of a qubit by a BEC produces a slow state- 
selective decoherence which obeys a power, rather that 
exponential, law in time (except for a static qubit, where 
the decoherence is extremely fast). Although this re- 
sult was obtained in the limit N — ► oo, it can be con- 
firmed by numerical evaluations of the qubit's density- 
matrix, Eq. © for a large but finite TV, Fig. It is 
this non-exponential relaxation and a strong dependence 
on the qubit's initial state that distinguishc the BEC 
model, with a single energy level in each of the reservoirs, 
from the exponential decoherence generated by a (Marko- 
vian) environment with a continuum spectrum of avail- 
able states. Common to both environments is, however, 
freezing of the qubit's internal transitions in the str ong 
interaction limit. This kind of Zeno effect [ljj, 120, [21| 
produced by the unitary evolution in the presence of an 
environment is somewhat different from its conventional 
prototype [22|, |23| which arises from frequent observa- 
tions of the evolving system. One remarkable feature of 
a Markovian environment is that the qubit's evolution 
under such frequent observations is practically indistin- 
guishable from its unitary observation- free evolution 24 1 . 
For a qubit-BEC hybrid system whose behavior is explic- 
itly non-Markovian, we expect the two types of evolution 
to be drastically different. A detailed investigation of this 
problem is, however, beyond the scope of the present pa- 
per. 
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